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D ! Abstract 



A study of W"'"W~ events accompanied by hard photon radiation produced in e^e^ colhsions 
at LEP is presented. Events consistent with two on-shell W-bosons and an isolated photon 
Vh ' are selected from 183 pb~^ of data recorded at i/i = 189 GeV. From these data, 17 W+W~7 
■ - candidates are selected with photon energy greater than 10 GeV, consistent with the Standard 
Model expectation. These events are used to measure the e"'"e~ -^ W"'"W~7 cross-section within 
a set of geometric and kinematic cuts, ^^^^wy = 136±37±8 fb, where the first error is statistical 
and the second systematic. The photon energy spectrum is used to set the first direct, albeit 
weak, limits on possible anomalous contributions to the W+W~77 and W+W~7Z° vertices: 
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where A represents the energy scale for new physics. 
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1 Introduction 

The W+W" cross-section has been precisely measured at LEP over a range of centre-of-mass 
energies [^ 0], and is well described by the Standard Model (SM) expectation 0. In this paper 
these measurements are extended to include also W"^W~ events with energetic photons in order 
to probe the modelling of electro-magnetic radiation in the W^W~ pair creation process. No 
previous measurements exist of three vector boson production when at least two bosons are 
massive. In addition, the W"*"W~7 final state may be sensitive to anomalous contributions to 
the SM W"*"W~77 and W"''W~Z°7 quartic gauge couplings, shown in Figure |I[ 




Figure 1: Standard Model W"*"W 77 and W"''W Z°7 quartic gauge couplings 



The non-Abelian nature of the electroweak sector of the Standard Model results in vector 
boson self-interactions. In addition to the triple gauge boson couplings (TGCs), W"'"W~7 
and W"'"W~Z°, the Standard Model predicts the existence of four quartic gauge couplings, 
W+W-W+W", W+W-Z°Z°, W+W-Z°7 and W+W-77. These couplings are not expected to 
play a significant role at LEP energies, but will be important at a future TeV linear collider Q] 
and at the LHC |^. There has been neither direct experimental confirmation of the existence 
of quartic couplings, nor any direct limits on possible anomalous quartic couplings. However, 



indirect limits on anomalous quartic couplings can be derived from the precise LEP/SLD Z° 
data [H, 0]. 

Quartic gauge boson couplings can be probed in final states with three vector bosons. At 
LEP centre-of-mass energies, final states involving three massive gauge bosons are kinematically 
out of reach. However, it is possible to study the processes e"'"e~ -^ W"'"W~7 and e"'"e~ —>■ Z°77. 
In the Standard Model, the contribution of the quartic coupling to e"'"e~ — * W'*'W~7 is expected 
to be too small to measure and that to e"'"e~ — * Z°77 is zero. Nevertheless, it is possible to set 
the first direct limits on possible anomalous contributions to the quartic gauge boson couplings. 

This paper describes a study of the process e+e~ — *> W+W~7 using 183 pb^^ of data 
recorded at a centre-of-mass energy of 189 GeV with the OPAL detector at LEP. It begins 
with an introduction to the phenomenology of quartic gauge boson couplings at LEP, followed 
by a description of the selection of W"'"W~7 events and of the measurement of the W+W~7 
cross-section for photon energies E^ > 10 GeV. The results are then interpreted as direct limits 
on possible anomalous W+W~77 and W+W^Z'^7 quartic gauge couplings and compared with 
the indirect limits. 

2 Anomalous Quartic Gauge Couplings 

In the Standard Model the form and strength of the vector boson self-interactions are fixed by 
SU(2)xU(l) gauge invariance. As is the case for triple gauge boson couplings P, in extensions 
to the Standard Model, anomalous quartic couplings can be parametrised by additional terms in 
the Lagrangian |^, TU, rT|. These are required to conserve custodial SU(2)c symmetry in order 
to avoid deviations of the p parameter^ from the experimentally well established value close to 
1. Only operators which do not introduce anomalous triple gauge couplings are considered. For 
example, the anomalous quadrupole moment operator generates both W+W~7 and W+W~77 
couplings. Therefore, it is not considered as a source of genuine anomalous quartic couplings 
since its strength, A-y, is already tightly constrained from the study of TGCs at LEP [|l], |1^ 
and at the Tevatron |13|. The lowest dimension operators which generate genuine anomalous 
quartic couplings involving photons are of dimension six. Three such possibilities are considered 
here, £0,£U| and £^1,111: 



with 






Z^/ cos9w 



where F^*^ and W^'^ are the field strength tensors of the photon and W fields respectively. Both 
£q and £g, which conserve C and P (separately), generate anomalous ^^^^^^"77 and Z''Z°77 

^p = Af-^/(M|o cos^ dw), where Mw and Afzo are the masses of the W^ and Z" bosons and 9w is the weak 
mixing angle. 



couplings. The CP violating term £g results in an anomalous W"''W~Z°7 coupling^. In each 
case, the strength of the coupling is proportional to a^/A^, where A represents a scale for new 
physics. A more general description of the operators leading to anomalous quartic couplings 



accessible at LEP can be found in the recent paper of Belanger et al. [|I4 . 

Figure ^ shows the predicted cross-section for W+W~7 as a function of different values of 
oq/A^, ttc/A^ and a„/A^, using the calculation of Stirling and Werthenbach |11|. The forms of 
/^g, Lq and £g dictate that contributions to the W+W~7 matrix element from the anomalous 



quartic couplings scale linearly with the energy of the photon ||TT| . The cross-section calculation 



of reference |TI| has been used to investigate further the effects of anomalous couplings on the 
energy spectrum of photons in W^W~7 events and on the angular distribution of the photons. 
For ^/s = 189 GeV the effects of a large anomalous coupling {e.g. a^) are shown in Figure ||b 
and ||c. The signal of anomalous quartic gauge boson couplings at LEP would be an excess of 
W"''W~7 events with photon energy greater than 10 GeV. 

3 The W+W-7 Final State 

There are 23 lowest order diagrams which contribute to the process e+e~ -^ W+W~(7) — *■ 
qqiVi'y. There are an additional 3 diagrams for the process e+e~ — > W+W^(7) -^ qqqq7- These 
diagrams correspond to initial state radiation (ISR), final state radiation (FSR), radiation from 
the W-boson (WR) and the Standard Model quartic gauge coupling (QGC) diagram. Different 
Monte Carlo generators and calculations employ different subsets/implementations of these 



diagrams. In this study, two Monte Carlo generators, WWF [T^ and KORALW IjT^ are used 
to simulate the Standard Model expectation for the W^W~7 final state. The WWF Monte 
Carlo generator produces W"'"W~7 events above some photon energy cut and within a specified 
acceptance in polar angle. This program implements a full matrix element calculation for the 
Standard Model ISR, FSR, WR and QGC diagrams [jl8|. However, this treatment may not be 



optimal for radiation from quarks where it is also necessary to consider QCD corrections to the 
radiation process. Higher order corrections are implemented using a structure function approach 
to coUinear ISR. The KORALW generator, used to simulate W+ W~ (7) events, does not include 
contributions from radiation from the intermediate W-bosons or from the Standard Model 
quartic coupling diagram. FSR from leptons is implemented using the PHOTOS [|19l program 



and FSR from quarks is modelled with JETSET ||20|. However, since the Standard Model 
W"'"W~7 cross-section is dominated by ISR, reasonable predictions are obtained. The Monte 
Carlo samples generated with both WWF and KORALW are passed through the full OPAL 



detector simulation E^. Neither WWF nor KORALW currently allows for anomalous quartic 



gauge couplings. The calculation of Stirling and Werthenbach [^ allows for the assessment 



of the impact of anomalous quartic couplings and is implemented in the EEWWG program. 
This calculation includes the ISR diagrams, the WR diagrams, the SM quartic diagram and 
can accommodate anomalous quartic couplings. It does not include FSR and therefore ignores 
any possible interference effects between the FSR diagrams and the other diagrams. 



^In reference ||lj| it is pointed out that the form for Cg used in [^ is missing a factor of i which is crucial 
for hermiticity. The missing factor of i was a typographical error in the paper. However, the implementation 
in the EEWWG program is correct ll^. 



3.1 W+W-7 Signal Definition 

The process e"'"e~ -^ W+W^7 results in a four-fermion plus photon final state, /1/2/3/47, 
where the fermion fiavours are appropriate for W-decay. Two distinct regions in phase space 
are considered. In the first, the /1/2 and fsf^ systems are produced with invariant masses 
close to Mw- In the second region, which is dominated by final state radiation, either the 
/1/27 and 73/4 combinations or the /1/2 and 73/47 combinations give invariant masses close 
to the W-boson mass. The second region has little sensitivity to anomalous quartic couplings. 
For photon energies E^ ^ Fw, where Fw is the W-boson width, it is possible to separate the 
phase space of events from FSR from those arising from the ISR or QGC diagrams using the 
invariant masses of the /1/2 and 73/4 systems. 

Both the selection procedure and the signal definition are designed to suppress contributions 
from FSR which is effectively a background process. The cuts used also reject photons from the 
decays of 7r° and rj mesons associated with the hadronic jets. For this paper, the W~*"W~7 — *■ 
/1/2/3/47 cross-section, denoted by 0"ww7) is defined by: 

• |cos^^| < 0.9, where cos 6'^ is the cosine of the polar angleQof the photon. 

• cosO^f < 0.9, where cos6-yf is the cosine of the minimum angle between the photon and 
any of the charged fermions in the four-fermion final state. 

• minimum(Mr y , Mf^j ) > 73 GeV, where Mf_-j: and Mf,-r are the invariant masses of 
the fermions from the W~ and from the W"*". 

• -E^ > 10 GeV, where E^ is the photon energy. 

In the above definition, all cuts are made on generator level quantities. Generator level refers 
to the true four-momenta of particles in the /1/2/3/47 final state. The first two cuts suppress 
ISR and FSR respectively and the invariant mass cuts further suppress final state radiation. 
The cross-section within the above kinematic cuts, (Tww7; is dominated by doubly-resonant 
W"'"W~ production. At ^/s = 189 GeV, the effects of other four-fermion diagrams resulting 
in final states which can interfere with W^W~7 have been evaluated with the RACOONWW 



program p2[. For the signal acceptance, the cross-section for the e"'"e — > W^W 7 diagrams 



alone is 109.33 ± 0.02 fb and when all (interfering) four-fermion diagrams are included this 



becomes 109.84 ± 0.03 fb |23], where in both cases the errors are statistical. Currently the 
predictions of RACOONWW do not include higher-order corrections from collinear ISR. 

The Standard Model expectations for the cross-section for the W"'"W~7 final state within 
this acceptance for E^ > 10 GeV have been determined using the WWF and KORALW gener- 
ators and the EEWWG program. The results are summarised in Table |I|. Also given in Table ^ 
are the cross-sections from EEWWG, WWF and RACOONWW where the effect of collinear 
ISR has been neglected. 

4 The OPAL Detector, Data and Monte Carlo Samples 

The OPAL detector includes a 3.7 m diameter tracking volume within a 0.435 T axial magnetic 
field. The tracking detectors consist of a silicon micro-vertex detector, a high precision gas 



•^The OPAL right-handed coordinate system is defined such that the origin is at the centre of the detector 
and the z axis points along the direction of the e^ beam, 9 is the polar angle with respect to the z axis. 



vertex detector and a large volume gas jet chamber. Lying outside the solenoid, the electro- 
magnetic calorimeter consisting of 11 704 lead glass blocks has full acceptance in the range 
|cos^| < 0.98 and a relative energy resolution of approximately 6% for 10 GeV photons. 
The magnet return yoke is instrumented with streamer tubes which serve as the hadronic 
calorimeter. Muon chambers outside the hadronic calorimeter provide muon identification in 



the range | cos9\ < 0.98. A detailed description of the OPAL detector can be found in |£4 . 

The integrated luminosity of the data sample recorded during the 1998 LEP operation, 
evaluated using small angle Bhabha scattering events observed in the silicon tungsten forward 



calorimeter pH], is 183.05±0.16(stat.)±0.37(sys.) pb~ . The luminosity-weighted mean centre- 
of-mass energy for the data sample is a/s = 188.63 ± 0.04 GeV. 

Both the WWF and KORALW Monte Carlo samples are used to simulate W^W~7 events. 
Selection efficiencies are estimated using WWF. Differences between efficiencies calculated 
using WWF and KORALW are assigned as systematic uncertainties related to the theoreti- 
cal modelling of W"*"W~7 production, and in particular, the modelling of FSR. A number of 
Monte Carlo samples, all including a full simulation of the OPAL detector, are used to simu- 
late background processes to W"^W~(7). The background process Z°— >qq is simulated using 



PYTHIA dg, with HERWIG ^ and KK2F p8| used as alternatives to assess systematic 



uncertainties. The four-fermion backgrounds from {Z^/'~f){Z^/'~f) and WeUe are simulated with 



PYTHIA, GRC4F |2| and EXCALIBUR |0|. 



5 W+W 7 Event Selection 

The selection of W'*"W^7 events proceeds in three stages: selection of W"''W~ events, photon 
identification and acceptance cuts, and kinematic requirements. Initially the energy require- 
ment on identified photons is E'^ > 5 GeV. However, for the determination of the W"'"W~7 cross- 
section and the anomalous coupling analysis only events with photon energies E^ > 10 GeV are 
used. The photon energy range 5-10 GeV has little sensitivity to anomalous couphngs and is 
subject to a larger background contribution and consequently larger systematic uncertainties. 
For this analysis, W^W~(7) — > tP^tPp'-y events are not used because it is not possible to 
reject adequately FSR photons using the kinematics of the event. However, in the determination 
of (Tww7 the value quoted has been corrected to include all final states. 

5.1 W+W- Selection 

The W"'"W~^qq£z7£ and W"^W~-^qqqq selections of reference [^ are used as the basis of the 
W"'"W~7 selection at 189 GeV. For the selection of W"'"W~(7) — >qq£z7£7 events, the standard 
W"'"W~^qq£l7£ selection is applied, but with one additional requirement. In a significant frac- 
tion of events selected as W^W^ -^ qq(T -^ h^mr^u)!/ and W'^'W^ -^ <1(i(t —>■ Sh'^mr^u)]/, the 
track(s) and clusters identified as the tau decay products are in reality associated with one of 
the jets. In order to reduce this contamination, the tau decay candidate track (highest momen- 
tum track in the case of 3-prong decays) is required to have momentum greater than 3 GeV. 
For photons with |cos^^| < 0.9 the selection efficiency for W"'"W~(7) ^qqiuij events is ap- 
proximately constant at 87% in the photon energy range 10-30 GeV, and is almost independent 
of the polar angle of the photon within the accepted region. 

For the selection of W"^W~(7) -^qqqqj events, a modified version of the W"''W~— >qqqq 
selection of reference |jl|] is used. In this selection events are forced into four jets using the 
Durham kr algorithm |^|. In approximately 10% of Monte Carlo events with a photon with 



E^ > 10 GeV, the photon alone forms one of the four jets. This introduces an additional 
inefficiency, due to the requirement in the preselection that there should be at least one charged 
track associated with each jet. For this reason, events failing the standard W+W~— s>qqqq 
selection are forced into four jets after excluding the highest energy isolated electro-magnetic 
calorimeter cluster and the selection re-applied. The selection efficiency for W^W~^qqqq7 
events with photon energies in the range 10-30 GeV is 90% and is approximately independent 
of the photon energy and photon polar angle. 

5.2 Photon Identification and Acceptance 

Photon identification is similar to that described in |32|. Photon candidates are identified as 
one of three types: 

• Unassociated electro-magnetic calorimeter clusters defined by the requirement that no 
charged track lie within the angular resolution of the cluster when extrapolated to the 
calorimeter. The lateral spread of the cluster was required to satisfy the criteria described 
in reference p2| . 



• Two-track photon conversions which are selected using an artificial neural network. 

• Conversions where only a single track is reconstructed, identified as an electro-magnetic 
calorimeter cluster associated with a charged track which is consistent with pointing to 
the primary vertex. The track is required to have no associated hits in either layer of the 
silicon micro- vertex detector or in the first 6 layers of the central vertex chamber. Up to 
one additional charged track passing the same criteria is allowed to point to the cluster. 

For both types of conversions, the photon energy is defined by the sum of cluster energies 
pointed to by the track(s). 

Photon candidates identified using the above criteria are required to satisfy isolation re- 
quirements. The energy of additional tracks and clusters in a 20° half-angle cone defined by 
the photon direction has to be less than 2 GeV. In addition, the energy deposited in the hadron 
calorimeter in a 20° half-angle cone around the photon candidate is required to be less than 
5 GeV. 

The identified photon is required to lie within the polar acceptance, 

• |cos^^| < 0.9. 

The photon is also required to be isolated from the charged fermions in the final state. Cuts 
are applied on the cosine of the angle between the photon and closest jet, cos 6*^- jet, and on 
the cosine of the angle between the photon and lepton, cos^^^: 

• cos 6*^- JET < 0.9, 

• cos^^^ < 0.9 for W"'"W^^qqez7e7 and W+W^^qqyuz7^7, 

• cosO^i < 0.7 for W+W"^qqrI7^7. 

For photons within the generator level acceptance E^ > 10 GeV, |cos^^| < 0.9 and cos6'^/ < 
0.9, the photon identification efficiency is about 83% for selected W"''W~(7) —>-qqiT'£'y and 
W"''W~(7) ^qqqq7 events. 
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5.3 Kinematic Requirements 

Kinematic cuts which reject events with FSR photons or with photons associated with the 
hadronic jets, i.e. from vr" and rj decays, are performed by comparing the results of three 
kinematic fits. In each case, the constraints of energy and momentum conservation are imposed. 
The finite W width is neglected and the two reconstructed masses of the W boson candidates 
are required to be equal |3^. The fits employed correspond to the following hypotheses: 



a) FSR(quark), assuming a two-body W^W final state, where the identified photon is 
associated with one of the jets. 

b) FSR(lepton), only used for W~''W~(7) — >qq£z7^7 events, assuming a two-body W+W~ 
final state, where the photon is associated with the charged lepton. 

c) ISR/QGC, assuming a three body final state consisting of two W bosons and the photon. 

An event is considered consistent with one of the above fit hypotheses if the fit converges with 
a fit probability of greater than 0.1% and if the reconstructed W boson mass is greater than 
74 GeV. Only W~^W~7 candidate events consistent with being two approximately on-shell 
W-bosons (fit c), M.j ~ M. j ~ Afw, are retained for the measurement of the W"^W~7 
cross-section and quartic gauge coupling study. This procedure suppresses events with final 
state radiation and events where the photon is from the decay of mesons. By reducing con- 
tamination from FSR, the sensitivity to anomalous quartic couplings is improved. In addition, 
the systematic uncertainties from photons associated with jets (FSR and vr^/r] decays) are 
significantly reduced. 

Selected W"'"W~(7) — ^qq£z7£7 events are required to be consistent with the ISR/QGC hy- 
pothesis using the above criteria. The effect of the invariant mass cut is indicated in Figure |^. 
In order further to reduce contributions from FSR, if the event is also consistent with the hy- 
pothesis of FSR from the lepton, it is required that the reconstructed mass from the ISR/QGC 
fit be closer to Mw than the mass obtained from the FSR(lepton) fit. 

In fully hadronic events there are three possible jet-pairing combinations. Only events where 
one of the three combinations is kinematically consistent with the ISR/QGC hypothesis are 
retained. If more than one combination satisfies this requirement, the fit yielding the highest 
probability is used. Events which are also consistent with the FSR(quark) hypothesis are 
rejected if the FSR fit has higher probability than the ISR/QGC fit and if MpsR < 86 GeV and 
-^iSR < 78 GeV, where Mfsr and Mjsr are the respective reconstructed masses. 

The application of the above kinematic requirements retains approximately 80% of W^W~7 
events within the signal definition of Section O whilst rejecting 85%-98% (increasing with the 



photon energy) of events with photons from FSR or from the decays of mesons. 

6 W+W~7 cross-section 

Using the selection criteria defined in the previous section, 17 W"'"W~7 events with E^ > 
10 GeV are selected compared to the Monte Carlo expectation^ of 13.2 events. Figure Ija 
shows the photon energy spectrum for the selected W"'"W~7 event sample at ^/s = 189 GeV. 
Also shown are the data for photon energies in the range 5 — 10 GeV. One event is observed 



^Here KORALW was used for the expectation from W+W (7) and PYTHIA was used to simulate the 
backgrounds from qq, Z^Z^and WeFc 



with E^ > 30 GeV, compared to the expectation of 0.1. This event occurs in the electro- 
magnetic calorimeter barrel/endcap overlap region where the energy resolution is relatively 
poor. Figure |b shows the distribution of cos 9^ for the events with E^ > 10 GeV and Figure |c 
shows the distribution of the cosine of the angle between the photon and nearest charged 
fermion in the event. Good agreement between data and Monte Carlo is observed. 

New physics could appear as resonant structure in the W7 invariant mass distribution. 
To investigate this possibility, the invariant masses of the two W^7 combinations in selected 
W"'"W~7 events are obtained from an additional kinematic fit. The fit uses the constraints of 
energy and momentum conservation and the constraint that the invariant masses of the recon- 
structed, /1/2 and 73/4 systems are both equal to the W-mass. Only events for which the kine- 
matic fit converges are retained. The W7 invariant mass is calculated from the four-momenta 
of the four fermions and the photon returned by the fit. Figure |^ shows the reconstructed 
invariant mass distribution for the two W^7 combinations for selected W"'"W~7 events with 
E^ > 5 GeV. No resonant structure is observed. 

The W"'"W~7 cross-section for E^ > 10 GeV is determined within the acceptance defined 
in Section p7T| . The selection efficiency for events within this acceptance, eww^, is evaluated 



using the WWF Monte Carlo sample. Backgrounds from qq(7) and Z^Z" are estimated using 
PYTHIA. The background from W+W" events with a fake photon from the electro-magnetic 
decays of mesons is estimated using KORALW. These three sources are summed to form the 
total non-W"'"W~7 background, ctbgd- Due to the steeply falling photon energy spectrum, 
the effect of the finite photon energy resolution is that more events from lower energies are 
reconstructed with E^ > 10 GeV than vice versa. A correction, Cres, is applied to account for 
this migration. In addition, it is necessary to account for accepted W''"W~7 events from outside 
the invariant mass region used to define the cross-section using a correction factor, Pww-y, 
which is defined as the fraction of accepted /i 72/3/47 events which have Mj:_j > 73 GeV and 
Mj:rF > 73 GeV. The W"'"W~7 cross-section is calculated from 

PwW-ri^ohs — O'BGD'^) 

(TWW7 — 7; , 

where A^obs is the accepted number of events with E^ > 10 GeV and £ is the integrated 
luminosity. The values of the quantities used to determine the cross-section and associated 
uncertainties are listed in Table 0, giving the result, 

aww7 = 136 ± 37 ± 8 fb, 

where the first error is statistical and the second systematic. This result is consistent with the 
expectations (85 - 102 fb) of the WWF, KORALW and EEWWG presented in Table |I|. The 
sources of systematic uncertainty, summarised in Table H, are described below. 

6.1 Systematic Uncertainties 

W+W" Selection Efficiency: The systematic uncertainty on the selection of W"'"W~ events 
is estimated to be 0.9%. This estimate was obtained in a manner similar to that described in 
[]I|. The largest uncertainties are related to the QCD and fragmentation modelling of jets. 

Photon Identification : A systematic uncertainty of 1% is assigned to cover the uncertainties 
in the simulation of the photon conversion rate and the accuracy of the simulation of the 
electromagnetic cluster shape. 
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Photon Isolation : The systematic error associated with the isolation requirements depends on 
the accuracy of the Monte Carlo simulation of the fragmentation process in hadronic jets. This 
is verified in T? — >■ qq events recorded at -^s = 91 GeV during the 1998 run. For each selected 
event, the inefficiency of the isolation requirements is determined for random orientations of 
the isolation cone and parametrised as a function of the angle between the cone and the nearest 
jet. For all jet-cone angles the inefficiency in the Monte Carlo and data agree to better than 
1%, consequently a 1% systematic error is assigned. 

Kinematic Fits : The W+W~7 event selections require that a kinematic fit converges and 
has a reasonable probability. Possible mis-modelling of the detector response/resolution could 
result in a difference in the rates at which the fits fail for data and Monte Carlo. This was 
checked by applying the kinematic fits used in the W mass analysis to all selected W"'"W~ 
events and comparing the failure rates for data and Monte Carlo. The difference in efficiency, 
1.2 ± 0.8%, is used as an estimate the systematic error associated with the kinematic fits. 

ECAL energy scale and resolution: The systematic error (2.5%) on the correction factor, 
Cres, which accouuts for feed-through from lower energy photons, is estimated by varying the 
electro-magnetic calorimeter energy resolution by ±25% and by varying the electro-magnetic 
calorimeter scale by ±1%. These variations were obtained by studying photons in e"'"e~ — * 77 
and e"'"e~ -^ e"'"e~7 at -\/i ~ Mgo and -/s = 189 GeV, and by studying the energy response to 
electrons. 

Modelling of photons from jets: The modelling of photon candidates associated with the 
hadronic jets (both from FSR and from vr*^ and rj decays) is studied by comparing the rate 
at which photons are identified in Z° —>■ qq events from the 1998 calibration data (y^=91.2 
GeV) to the PYTHIA prediction. For 10 GeV< E^ < 20 GeV, there are 40% more photon 
candidates identified in the data than expected from the Monte Carlo. Above 20 GeV the data 
are consistent with the Monte Carlo expectation. The ratio of data to Monte Carlo is used to 
estimate an energy- dependent correction to the Monte Carlo expectation for the background 
from W"'"W~ events with fake photons. The size of this correction is assigned as a systematic 
uncertainty. 

ISR Modelling in qq7 Background: The dominant source of non-W^W^ background is 
from e"'"e^ -^ 2^/7 -^ qq7 where the identified photon is a genuine photon from ISR. The 
accuracy of the simulation of initial state radiation with |cos6'^| < 0.9 is studied using identified 
ISR photons in multi-hadronic events [Q recorded at y/s = 189 GeV. A kinematic fit is used 



to reject events with additional radiation along the beam direction. In the range of photon 
energies of relevance to the W'^W~7 analysis, 5 GeV < E^ < 50 GeV, the fractions of events 
with an identified photon in data and Monte Carlo agree to better than 20% which is assigned 
as a systematic uncertainty. 

Four-fermion Background: For the W^W~7 cross-section measurement PYTHIA is used 
to estimate the background from Z°Z'^ and single W production (Wez/g)- This four-fermion 
background estimate was compared to that obtained using GRC4F and EXCALIBUR. The 
largest difference is used to assign a systematic error. 

W+W~7 events from outside the signal definition: In order to account for accepted 
W^W^7 events from outside the cross-section signal definition, a correction factor is used. 
This factor, Pww7, is determined to be (85.3 ± 0.7)% from WWF and (85.0 ± 1.0)% from 
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KORALW. The WWF value is taken as the central value. The comparison of these generators 
is sensitive to uncertainties in the treatment of FSR. Although no discrepancy is observed, the 
statistical error on the difference is taken as a systematic error (1.2%). In addition, a systematic 
uncertainty of 1.5% is assigned to account for the differences between data and Monte Carlo in 
the modelling of photons from jets, described previously. 

7 Limits on Anomalous Quartic Couplings 

The selected W"'"W~7 events are used to set limits on possible anomalous contributions to the 
W"'"W~77 and W"'"W~Z°7 quartic gauge couplings. The limits are extracted from the measured 
differential cross-section as a function of the photon energy and photon polar angle. 

There is currently no Monte Carlo available that implements anomalous quartic couplings 
and has contributions from all Standard Model diagrams. The expectated differential W^W~7 
cross-sections are obtained using the EEWWG program which allows for anomalous quartic 
couplings. Since FSR is not implemented in EEWWG it is included as an additional background 
process. The Monte Carlo treatment of signal and background is described below. 

7.1 W+W-7 Signal 

The EEWWG calculation is used to determine the expected contribution from the ISR, WR, 



Standard Model QGC and anomalous QGC diagrams within the signal definition of Section [O. 
The only detector effect that is included is a Gaussian smearing associated with the energy 
resolution for the photon. No smearing was applied to the photon or fermion angles since 
these are relatively well measured. The selection efficiency is estimated (for events within 
the signal acceptance cuts) from the WWF Monte Carlo and is parametrised as a function of 
photon energy. The small dependence of the selection efficiency on cos 9^ is neglected in the 
determination of the limits on anomalous couplings. The validity of this procedure is verified 
using a small sample of fully-simulated Monte Carlo events generated using the EEWWG 
calculation. 

The contribution to the cross-section from any anomalous coupling increases rapidly with 
centre-of-mass energy due to the additional phase space for higher energy photons. Thus the 
contribution to the W"'"W~7 cross-section from anomalous QGCs is sensitive to higher-order 
QED radiative corrections. The EEWWG calculation currently only includes W~''W~7 final 
states. For the case of the QGC diagram, it therefore neglects the effect of additional ISR 
which reduces the average effective centre-of-mass energy, ys' . The effect of collinear ISR from 
both the electron and positron has been incorporated into the EEWWG calculation using the 
collinear radiator function from EXCALIBUR. At ^/s = 189 GeV, QED radiative corrections 
reduce the cross-section from anomalous QGCs by 26%. 

7.2 Background Treatment 

The contribution from FSR is estimated using KORALW, with FSR photons being identified 
using the Monte Carlo generator information. KORALW is also used to estimate the contri- 



bution from other diagrams producing final states outside the signal definition of Section |3?1 . 
Neglecting the contribution of anomalous couplings to events outside the signal definition not 
only results in conservative limits but reduces the possible effects of interference with Stan- 
dard Model FSR diagrams. Background from W"''W~ with a fake photon is calculated using 
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the corrected background from KORALW. Similarly, the four-fermion and qq backgrounds 
are estimated using PYTHIA. Table ^ summarises the background contributions used in the 
anomalous coupling study. In the table, the W"'"W~7 background has been divided into that 
from ISR (outside the signal definition mass and/or acceptance cuts) and that from FSR. 

7.3 Results 

For the case of Standard Model quartic gauge couplings, the expected photon energy spectrum 
for selected events estimated using the above procedure is shown in Figure ^. The spectrum 
is in reasonable agreement with that obtained from the WWF generator. By combining the 
predicted cross-sections from EEWWG and the backgrounds from Table ^ the expected distri- 
butions corresponding to non-zero anomalous coupling parameters can be generated. Figure ^ 
shows the expected photon energy distributions for two values of the anomalous coupling oq/A^ 
(in addition to the Standard Model prediction). 

To set limits on possible anomalous couplings a binned maximum likelihood fit to the 
observed distribution of [E^, cos 6'^] is performed using bins of [5 GeV, 0.1], with the fit range 
restricted to E-^ > 10 GeV. The 95% confidence level upper limits on the anomalous couplings 
are obtained from the resulting likelihood curves: 

-0.070 GeV-^ < ao/A^ < 0.070 GeV"^ 
-0.13 GeV-^ < ajk^ < 0.19 GeV-^ 
-0.61 GeV"^ < a„/A2 < 0.57 GeV"^ 

The expected 95% confidence level upper limit on ao/A^ is 0.045. These limits from the data 
are higher than the expected limits due to the slight excess of high energy photons in the data. 
The probability of obtaining a limit greater than or equal to the observed limit is approximately 
2%. The limits were derived including the systematic uncertainties described in the previous 
section. For these results, a 15% theoretical uncertainty on the cross-section normalisation of 
the EEWWG calculation is assumed. 

These are the first direct limits on anomalous quartic couplings. However, indirect limits 
can be derived from radiative corrections to LEP/SLD Z^ data. The experimental limits on 
deviations of the S,U,T parameters [^ (or equivalently ei, €2, £3 [Q) from their Standard Model 
values has been used to place constraints on the anomalous W^W^77 and Z°Z°77 couplings 0: 
-0.0007 < ao/A^ < 0.0001 and -0.0017 < aJA^ < 0.0009. However, these indirect limits are 
obtained under a restrictive set of assumptions and their validity has been questioned fl^ . 
There are currently no limits on a„ corresponding to an anomalous W^W~Z''7 coupling. 

8 Conclusions 

The first study of W"'"W~7 events produced in e"'"e~ collisions is presented. From the seventeen 
selected events with E^ > 10 GeV, the W"'"W~7 production cross-section is measured to be: 

aww7 = 136 ± 37 ± 8 fb, 

within the /1/2/3/47 generator level acceptance defined by |cos6'^| < 0.9, cosO^f < 0.9 and 
Mr J , Mr J > 73 GeV, in agreement with Standard Model expectation. 
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The distribution of the photon energy and polar angle is used to set limits on possible 
anomalous contributions to the W"'"W~77 and W"'"W~7Z° vertices: 



-0.070 GeV" 


2< 


ao/A^ 


< 


0.070 GeV-^ 


-0.13 GeV- 


-2< 


flc/A^ 


< 


0.19 GeV-2, 


-0.61 GeV" 


-2< 


a„/A2 


< 


0.57 GeV-^ 



where A represents the energy scale for new physics. These are the first direct limits on 
anomalous quartic couplings. 
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Program 


Collinear ISR 


^WW7 


KORALW 

EEWWG 
WWF 


yes 
yes 
yes 


102.4 ±2.2 fb 
85.5 ± 0.4 fb 
89.3 ± 0.6 fb 


EEWWG 

WWF 

RACOONWW 


no 
no 
no 


106.1 ± 0.4 fb 

110.3 ± 0.6 fb 

109.33 ±0.02 fb 



Table 1: Standard Model cross-sections for the process e+e~ -^ W"'"W~7 within the /1/2/3/47 
generator level acceptance defined by E-y > 10 GeV, |cos6'-y| < 0.9, cos 9^f < 0.9 and 
Mr -J , Mr ^ > 73 GeV. Where possible, values are quoted corresponding to the cases with 



'/1/2 



hf 



and without collinear ISR. 



Quantity 


Value 


A^obs 


17 


C 


183.1 ± 0.6 pb-i 


^WW7 


48.7 ±1.8% 


Crcs 


1.065 ±0.026 


0"BGD 


9.9 ±2.0 fb 


PwW7 


85.3 ±1.9% 



Table 2: Values of the quantities used in the determination of the W"'"W 7 cross-section. The 
errors include components from Monte Carlo statistics and systematic uncertainties. 
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Source 


A0-WW7 


W+W" Event Selection 


3.0 fb 


qq7 efficiency 


0.5 fb 


Photon Identification 


1.6 fb 


Photon Isolation 


1.6 fb 


Kinematic Fits 


1.6 fb 


ECAL energy scale 


1.8 fb 


ECAL energy resolution 


2.8 fb 


Photons from jets 


2.7 fb 


ISR modelling qq7 


1.9 fb 


Acceptance Cuts 


2.8 fb 


Four-fermion Events 


2.0 fb 


PwW7 


3.0 fb 


Total 


7.7 fb 



Table 3: Contributions to the systematic error on awwi- 



-E-y Range 


non-W"''W 7 Background 


W+W 7 Background 




W+W" 


4/ 


qq 


ISR 


FSR(lepton) 


FSR(quark) 


10 - 15 GeV 


0.36 ±0.08 


0.08 ±0.03 


0.32 ±0.08 


0.29 ±0.05 


1.66 ±0.30 


0.74 ±0.14 


15 - 20 GeV 


0.19 ±0.06 


0.10 ±0.04 


0.26 ±0.07 


0.25 ±0.05 


0.15 ±0.04 


0.36 ±0.06 


20 - 25 GeV 


0.05 ±0.03 


0.02 ±0.01 


0.16 ±0.05 


0.06 ±0.03 


0.01 ±0.01 


0.01 ±0.01 


25 - 30 GeV 


0.04 ±0.02 


0.04 ±0.02 


0.07 ±0.03 


0.00 ±0.01 


— 


0.11 ±0.04 


30 - 35 GeV 


0.00 ±0.01 


— 


0.02 ±0.02 


— 


— 


0.01 ±0.01 



Table 4: Expected numbers of background events used in the anomalous coupling analysis. 
The W"'"W~ background consists of events with a photon from the decays of a tt*^, t] etc. In 
addition, W"'"W~7 events from FSR diagrams are considered as background and estimated 
using KORALW. Accepted W"'^W~7 events from outside the invariant mass cuts of the cross- 
section signal definition are also treated as background (denoted ISR). The quoted errors include 
systematic contributions. 
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Figure 2: The effect of anomalous quartic couplings in e^e~ collisions at ^/s =189 GeV. a) 
Expected W"'"W~7 cross-section for different anomalous quartic couplings, the photon was 
required to have energy, Ey > 10 GeV and be within the polar acceptance, |cos^-y| < 0.9. 
The kinematic cuts used to define the cross-section, see text, were also applied. Plots b) and 
c) show the effect of an anomalous quartic gauge coupling, ao, on the energy and angular 
distribution of photons in W^W"7 events, where the fiducial cuts, E^ > 5 GeV and |cos^^| < 
0.9 were imposed. In all three figures, the contribution from final state radiation diagrams is 
not included. 
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Figure 3: The reconstructed W-boson mass distribution for candidate W"'"W~(7) — s>qq£z7^7 
events using the ISR/QGC hypothesis. The points show the data. The histogram shows the 
Monte Carlo expectation where KORALW has been used to simulate the W"'"W~ (7) final state 
and PYTHIA is used to simulate non-W"'"W~ background. The hatched histogram shows the 
contribution from W"'"W~ events with final state radiation or where the photon candidate 
originates from the decay of a meson. The double-hatched histogram shows the contribution 
from non-W"'"W~ background. The cut is indicated by the arrow. 
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Figure 4: For selected W"'"W~7 events, a) shows the photon energy spectrum, b) the cosine 
of the polar angle of the photon, and c) the cosine of the angle between the photon and the 
nearest fermion. In b) and c) the cut E^ > 10 GeV has been applied. The 189 GeV data are 
shown by the points with error bars and the total SM Monte Carlo expectation shown by the 
histogram. The hatched histograms indicate the contribution from non-W"'"W~ background 
and from fake photons in selected W"'"W~ events. In Figure a) the predictions from KORALW 
are compared to those from WWF and EEWWG. The prediction of EEWWG also includes 
the estimated contribution from FSR obtained using the KORALW generator. 
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Figure 5: Reconstructed invariant mass of W^7 in selected W^W~7 events with E^ > 5 GeV 
(two entries per event). The 189 GeV data are shown by the points, the Standard Model 
expectation, determined from KORALW and PYTHIA is shown by the histogram. The shaded 
histogram shows the non-W"'"W~ (7) background. Also shown, by the dotted histogram, is the 
expected reconstructed mass distribution (arbitrary normalisation) for a narrow resonance of 
mass 100 GeV which decays to W7. 
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Figure 6: Energy spectrum of photons in candidate W"'"W~7 events. The data, shown by the 
points with error bars, are compared to the expected distributions for three different values 
of the anomalous coupling a^/h?; the Standard Model value a^/ h? = 0, the best fit value 
ao/A^ = 0.044 and the distribution corresponding to the 95% confidence limit value of ao/A^ = 
0.070. The double hatched histogram shows the background from non-W^W~ events and from 
W"'"W~ events with a fake photon. The singly hatched histogram shows the contribution from 
final state radiation. 
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